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ABSTRACT 

Using the stellar photometry catalogue based on the latest data release (DR4) of the 
Sloan Digital Sky Survey (SDSS), a study of the Galactic structure using star counts 
is carried out for selected areas of the sky. The sample areas are selected along a 
circle at a Galactic latitude of +60°, and 10 strips of high Galactic latitude along 
different longitudes. Direct statistics of the data show that the surface densities of I 
from 180° to 360° are systematically higher than those of I from 0° to 180°, defining 
a region of overdensity (in the direction of Virgo) and another one of underdensity (in 
the direction of Ursa Major) with respect to an axisymmetric model. It is shown by 
comparing the results from star counts in the {g — r) colour that the density deviations 
are due to an asymmetry of the stellar density in the halo. Theoretical models for the 
surface density profile are built and star counts are performed using a triaxial halo of 
which the parameters are constrained by observational data. Two possible reasons for 
the asymmetric structure are discussed. 
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1 INTRODUCTION 

It has been one of the central goals of human knowledge to 
understand the position of the Earth in the universe. A lot of 
effort has been spent to derive the structure and the history 
of formation and evolution of our Galaxy. Doing this is a 
particularly difficult task, simply because there is no good 
viewpoint, as opposed to the case for external galaxies whose 
structures can be observed directly. The structure, formation 
and evolution of the Galaxy are very important issues in 
contemporary astrophysics; they are closely related to stellar 
quantities such as distance, age, metallicity, and kinematics 
of the individual stars, among others. These properties can 
be obtained for large stellar samples using straightforward 
photometric and spectroscopic observations. In the present 
era, before the full exploitation of the planned huge spectral 
surveys of Galactic stellar objects (e.g. GAIA, SEGUE, 6dF, 
LAMOST) is possible, using star counts based on all-sky 
photometric surveys is a direct and one of the few accessible 
methods for the study of the structure of the Galaxy. 

Probing the Galaxy using star counts has a history of 
about 200 years (Reid 1993), with the modern epoch marked 
by the classic work of Bahcall & Soneira (1980), based on 
several observations (Seares et al. 1925, Kron 1978, Tyson 
& Jarvis 1979, Peterson, Ellis & Kibblewhite 1979). Bahcall 
& Soneira built a standard model, in which the Galaxy was 
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simplified and parameterized by an exponential disk and a 
spheroid characterized by a de Vaucouleurs profile. Assum- 
ing a density profile for each population, based on detailed 
studies of nearby spiral galaxies, and assuming a luminos- 
ity function as observed for stars in the local volume, their 
model predicts number density - magnitude - colour distri- 
butions. 

Further studies on this subject provided certain modifi- 
cations to the Bahcall & Soneira model. Theoretical models 
were further constrained using new observations were subse- 
quently obtained. In 1984, Gilmore claimed for the first time 
that there is a thick disk in our Galaxy. The star counts 
approach to the study of the structure of the Galaxy was 
extensively reviewed by Reid (1993); he, and his colleagues, 
(Reid 1993, Reid, Yan & Majewski 1996, Majewski, Siegel 
& Kunkel et al. 1999, Siegel, Majewski & Reid 2002) car- 
ried out a series of studies appropriately termed "star counts 
redivivus", dedicated respectively to: faint magnitude star 
counts, the halo luminosity function, and an exploration of 
the contamination of star counts by star streams. When the 
Sloan Digital Sky Survey (SDSS) made its first data release 
publicly available, it raised very promising prospects for the 
study of Galactic structure thanks to the large sky coverage 
and the deep stellar photometry. Chen et al. (2001) analysed 
the SDSS early data release (EDR); they confirmed the exis- 
tence of a Galactic thick disk and obtained a set of structural 
parameters for the Galaxy. With the increasing sky cover- 
age in subsequent SDSS data releases, it is expected that one 
will be able to determine the Galactic structure parameters 
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ever more precisely. Except for constraining the density pro- 
file, the stellar luminosity function can also be determined 
with a higher accuracy. In Gould, Bahcall & Flynnet (1996), 
and Zheng, Flynn & Gould (2001, 2004) the disk M-dwarf 
and 7-band luminosity functions are obtained using Hubble 
Space Telescope star counts. In Reid (2002), the luminosity 
function in the solar neighborhood was studied based on the 
Palomar/MSU nearby star spectroscopic survey. 

The Galactic structural parameters derived from vari- 
ous studies show large divergence, which cannot simply be 
attributed to statistical errors, although these are generally 
significant. In star counts from fields selected following a 
certain symmetry approach, there are deviations from the 
axisymmetric standard model which are obvious and cannot 
be neglected. Spergel & Blitz (1988) pointed out that this 
might be an indication of a triaxial structure of the Galaxy, 
and they used a triaxial halo to explain the asymmetries 
detected from gas motions around ^=180°. Later, Blitz & 
Spergel (1991) showed that the observed Hi gas is moving in 
an asymmetric gravitational potential. Using areas selected 
from the ASP catalogue, Larsen & Humphreys (1996) found 
deviations of star counts from the axisymmetric model; they 
found that there are significant excesses of stars in quad- 
rant I, at / = 20° — 50°, compared to quadrant IV, around 
b — 30° . Further work was published by Parker, Humphreys 
& Larsen (2003), who compare star counts in 120 POSS I 
fields. They claimed that the excess is not associated with 
any known stream such as those linked to the Sagittarius 
dwarf, Magellanic stream or Fornax-Leo-Sculptor stream. 
Parker, Humphreys & Beers (2004), however, concluded that 
the excess is due to thick disk stars and may be associated 
with the recently discovered Canis Major debris stream. In 
Hartwick (2000), possible reasons for causing triaxiality due 
to nearby spiral and dwarf galaxies are considered. Newberg 
& Yanny (2005b) built triaxial ellipsoidal stellar halo models 
and analysed F turn-off star samples from SDSS DR3. 

To explain the observed excess of the stellar number 
density with respect to an axisymmetric halo, two compet- 
ing scenarios have been proposed: the triaxialilty could be 
an inherent property related to the smooth structure of the 
Galaxy, or it could be related to the remnants of some his- 
torical merger events. Recently, this debate has intensified 
with the reported "Virgo overdensity" by Juric et al.(2005), 
which may be the cause of the large-scale asymmetry in the 
structure of the Galaxy. By using the photometric parallax 
method, they measured stellar distances which allowed them 
to construct a map of the 3D structure of the sky area cov- 
ered by current SDSS survey. Their study has, for the first 
time, revealed large-scale patterns in the Galaxy. In the same 
data that showed the Virgo overdensity, an underdense area 
with respect to an axisymmetric halo model is also obvi- 
ous in the direction of UMa. Whereas the overdensity can 
be naturally interpreted as a large-scale star stream in the 
halo model, the UMa underdensity cannot be accounted for 
in the same way. We believe that this observational fact is 
linked to the smooth structure of the halo in which the star 
streams (including the Virgo overdensity) are embedded. To 
tackle this problem, a similar star count analysis of Galactic 
structure will be followed in this paper. 

The Sloan Digital Sky Survey is a five-colour 
{u, g,r,i, z) broad-band photometric survey with a wave- 
length coverage from 3000 to 11,000A (Fukugita, Ichikawa 



& Gunn 1996). The DR4 imaging catalogue covers 6670 
square degrees. Its detection repeatability is complete at a 
95% level for point sources brighter than the limiting appar- 
ent magnitudes of 22.0, 22.2, 22.2, 21.3 and 20.5 for u,g,r,i 
and z, respectively (http://www.sdss.org). It has facilitated 
great progress in the study of the structure of the Galaxy, 
which was reviewed in the SDSS-II SEGUE Project. For the 
study of Galactic structure using star counts, SDSS provides 
the most up-to-date and complete stellar sample in terms of 
both sky coverage and magnitude limit. Thus, the SDSS is 
suitable to study the large-scale asymmetry of the Galaxy. 
In this paper, we build a model of the Galaxy with a non- 
axisymmetric geometry, and constrain the parameters of the 
model using observations in a set of characteristic directions 
taken from SDSS DR4. In Section 2, the selected observa- 
tional data are described. In Section 3, the model used to 
fit the observational data is introduced. The results of triax- 
ial halo model and the asymmetric structure of the Galaxy 
are described and analysed in Section 4. A discussion of the 
present results and our conclusions are given in the final 
section. 



2 OBSERVATIONAL DATA 

The SDSS is based on observations obtained with a 2.5m 
telescope with a 3° field of view at the Apache Point Ob- 
servatory in New Mexico, USA. Standard operations of the 
SDSS were performed from 2000 to 2005. Eventually, it wiU 
cover TT steradians centred on the northern Galactic cap and 
three stripes near the southern Galactic cap (Stoughton, 
Lupton & Bernardi 2002). The operating mode of observa- 
tions is by means of drift scanning. One continuous scan 
is designated by a "run", and two such runs generate a 
"stripe"; they have an overlap region of 1 arcmin on the 
edges of the two scans (Chen et al. 2001). Catalogues of ob- 
jects classified as "stars" from DR4 were downloaded from 
the sky server for sky areas selected to meet our aims. The 
retrieved data is free of overlap and has been deblended; 
saturated stars and poor-quality objects near bright stars 
and the edges of the frames were all filtered out (Newberg 
& Yanny 2005b). 

If the Galaxy is axisymmetric, as usually assumed in 
previous studies, the projected stellar number density at the 
same Galactic latitude will increase towards the direction of 
^ = 0°, and decrease when approaching £ = 180°, and any 
sky area pair mirror-symmetrically selected on either side 
of the ^ = 0° meridian plane should retrieve the same star 
number counts, within the uncertainties. If the axisymmet- 
ric model is correct, such a picture should be reproduced by 
star counts across the whole sky area, and can in principle 
also be detected in symmetrically selected subsets of sky ar- 
eas. Limited by the current SDSS sky coverage, whole-sky 
star counts are not possible, but a uniform data set covering 
all longitudes at a given latitude should suffice, as this would 
be a perfect probe to check the axisymmetric assumption. 
Such a data set is now available from SDSS DR4, which will 
be addressed in the following. Fig. 1 is a Lambert projection 
of the northern Galactic hemisphere, showing the selected 
sky areas used in this work (marked by diamonds). Each of 
the selected sky areas has a size of about 2° x 2° ; the Galac- 
tic coordinates of the centres of the selected fields are given 



The asymmetric structure of the Galactic halo 3 



in Table 1. The sky areas selected along the circle at 6 = 60° 
form the first group of our working data set. To check the 
structure at latitudes different from b = 60°, sky area pairs 
at different latitudes (from b = 55° — 85° , equally spaced by 
5°) along our longitudinal grid (in steps of 30° from to --^ 
360°) are also selected. In order to minimise the effects of 
extinction, sky areas below 6 = 55° were excluded. Prom the 
present SDSS public data release, sky areas of 11 longitudi- 
nal strips are available, and selected field from these form the 
remaining 11 data groups in our sample. From the second to 
the twelfth group, the sky areas within each group have the 
same longitude, so that the structural properties refiected 
by the stellar number densities along a given longitude can 
be derived. In addition, through comparison between sky 
areas having mirror symmetry with respect to the i = 0° 
meridian deviations from the axisymmetric Galaxy model 
can also be checked. Owing to the SDSS observing strategy, 
stars brighter than 15™ will be saturated, and star counts 
will be not complete for magnitudes fainter than 22™. 2 in 
the g and r bands, so that our work will be limited in the 
range of magnitudes from 15™ to 22™ in both the g and r 
bands. 



3 THE THEORETICAL MODEL 

To fully reveal the structure of the Milky Way, accurate 
measurements of the distances to individual objects are key 
ingredients. Such information is still lacking at the present 
time, however. As a consequence, the structure of the Galaxy 
cannot be resolved satisfactorily. Very large-scale photomet- 
ric surveys greatly improve our knowledge of Galactic struc- 
ture. However, the problem remains even if we use the data 
generated by the SDSS project, because of our inability to 
determine stellar distances accurately. The major science 
drivers of the SDSS arc focused on extragalactic objects. 
Therefore, a wide-filter photometric system was designed to 
detect the faintest possible extragalactic objects. The nar- 
rowest filter, u, is 600A. Such a photometric system has an 
obvious disadvantage when one wants to measure absolute 
stellar magnitudes, compared with the classical ubvyP sys- 
tem. The latter system by itself is already characterised by 
a scatter of 0.3 mag for A and F-type stars (Crawford 1975, 
1979). The difficulty to distinguish between giant and main- 
sequence stars will affect counts of G and K-type stars even 
more seriously in the SDSS system if one applies the colour- 
magnitude relation for main-sequence star to all stars. Fur- 
thermore, the large range of metallicity that is intrinsically 
present in the sample will inherently lead to large scatter 
in the colour-magnitude relation. To obtain distances to 
stars and to determine the structure of the Galaxy directly 
from SDSS data is thus somewhat difficult. A preliminary 
approach to this problem using the photometric parallax 
method and its pitfalls are carefully analysed by Juric et al. 
(2005). 

To derive the structure of the Galaxy from the inte- 
grated stellar light distribution along the line of sight, a 
theoretical model is needed. Our Galactic disk models fol- 
lows is based on that of Bahcall & Soncira (1980), while 
we have also included a thick-disk component. Because of 
large deviations from an axisymmetric structure (described 
in Section 4), we have modified the standard axisymmetric 



halo; these deviations will be accounted for by a triaxial halo 
instead. 

The basic star count equation we use is: 
A{Mi,M2,t,b) = / AM I R'^dRp{r)(P{M)dn, (1) 

J Ml Jo 

where A{Mi,M2,£, b) is the projected surface density in the 
absolute magnitude bin from Afi to Ad2 in the direction of 
Galactic longitude £ and Galactic latitude b; R is the helio- 
centric distance of a given star; p(r) is the density profile 
of each stellar population as a function of r, the distance to 
the Galactic Centre; </>(M) is the luminosity function of each 
contributing stellar population, which is a function of abso- 
lute magnitude and dQ is the solid angle element covered by 
the observations. 

For the thin and thick disk components, the density 
profile is assumed to be of the following exponential form, 

p(r)=exp[-|a:|/if-x//i], (2) 

where z is the height above the Galactic plane, x the pro- 
jected distance to the Galactic Centre in the Galactic plane, 
and H and h are the scale height and scale length in z and 
X, respectively. As this study focuses on the structure of the 
halo, we used fixed models for both the thin and the thick 
disks, so that complications due to extra free parameters for 
these disks can be excluded. The parameters of the disks are 
taken from Chen et al. (2001). For the thin disk, H = 330 pc 
and h = 2.25 kpc; for the thick disk, H = 650 pc and h = 3.5 
kpc. The galactocentric distance of the Sun is ro = 8.5 kpc, 
and the vertical distance of the Sun from the Galactic disk 
plane is zq = 27 pc. 

For the halo, a triaxial model is adopted to fit the ob- 
servational data. A set of suitable coordinate frames is built 
following Newberg & Yanny (2005b), of which the schematic 
diagram is shown in Fig. 2. Adopting the position of the Sun 
as the origin of the coordinate frame, the direction from the 
Sun to the Galactic Centre is taken as the x-axis. The y- 
axis is the direction perpendicular to the x-axis in an anti- 
clockwise sense within the Galactic plane, while the 2:-axis 
coincides with the direction towards the northern Galactic 
pole. The x' and y' directions are parallel to x and y, but 
their origin is at the Galactic Centre instead; these a^xes also 
define a right-handed coordinate system. Furthermore, R is 
the distance from a star to the Sun, r the distance from 
the star to the Galactic Centre and t is the projection of 
r on the Galactic plane. Therefore, z = R x sin(6) -I- zo, 
Zq = 27 pc, t = ^yrQ + R^ — 2 X Rx ro X cos(6) x cos(Z). 
The Galactic Centre coordinate frame is rotated in such a 
way that its axes are settled right on the three axes of the 
triaxial ellipsoid. The transformation between the cooordi- 
nate frames can be separated into three rotational steps: 
first, fix arbitrarily one of the three axes and rotate the 
other two. The rotation angles of the three steps, denoted 
by 0, (, and are defined when fixing the z, x and y 
axes, respectively. In Fig. 2, for example, where we rotate x' 
and y' around the z' axis, the relevant transformations are, 
Xi = x' X cos{d) + y' X sin(6l), yi = —x' x sm(d) + y' x cos{d), 
and zi = z' . The remaining two steps involve doing the 
same for the other two axes. The new axes obtained after 
the third rotation are (xa, j/3, 2:3), and the three axes of the 
triaxial ellipsoid are a, h and c, where o is the longest and c 
the shortest axis. The axis ratios of the triaxial ellipsoid are 
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p = b/a, and q = c/a, while r = + {ys/p)'^ + (z-i/q)'^. 
The axisymmetric halo is described hj p = 1. 

A power-law (Reid 1993) is adopted for the halo density 
distribution, 

PW = ^n^> (3) 

where ao = 1000 is a normalisation constant. Wc adopt the 
V-band luminosity function of Robin & Cre'e (1986) in our 
model, which provides sufficiently accurate two-dimensional 
distributions in both luminosity and spectral type down to 
V = 12 mag. A two-dimensional 1^-band luminosity func- 
tion can easily be transformed to that in the g band. Robin 
& Crezc (1986) compiled a number of luminosity functions 
available at the time, including those of Wielen (1974, 1983) 
and Houk & Cowley (1975, 1978, 1982). The work of Wie- 
len is widely used in star count studies to obtain number 
density-magnitude relations for 2 ^ A4v ^ 12 mag, while 
Houk & Cowley's luminosity function is often used for the 
bright stars (—2 ^ My ^ 2 mag). Deutschman, Davis 
& Schild (1976) and Hayes (1978) provided calibrations of 
stellar luminosity and spectral type. To deal with the data 
used for the present study, a further transformation relation 
(Chenet aL 2001) from V tog,g = "K-h0.53(B-y) -0.075, 
was also adopted. 

By using the density profile and the luminosity function 
discussed and presented above, the distribution in absolute 
g magnitude of the stellar projected surface number density 
can be obtained from Eq. (1). 

A colour-magnitude relation (CMR) is needed to obtain 
the number density distribution in the {g — r) colour. It is 
usual practice to use the CMR defined by colour-magnitude 
diagrams (CMDs) of star clusters. However, there are usu- 
ally too few stars in cluster CMDs with absolute magnitudes 
fainter than lO"*. Therefore, at fainter magnitudes we do not 
have a precise CMR. As an alternative approach, wc adopt 
the theoretical isoclirones of Girardi, Grebel & Odenkirchen 
(2004) to derive the CMR required here. The mean ages 
of thin disk, thick disk and halo are assumed to be 4.5, 11 
and 13 Gyr, and their metallicities are assumed to be, re- 
spectively, 0.019, 0.004 and 0.0012. The theoretical CMR is 
preferred because it includes faint stars that are otherwise 
not available from observations. However, we have to bear 
in mind that it cannot mimic the scatter on the relation in 
real situations. 

Given a theoretical model for stellar distributions of dif- 
ferent Galactic components as described above, it is straight- 
forward to calculate the projected surface number density 
in absolute magnitude bins in any given direction. From ob- 
servations, the information that can be directly extracted is 
the projected surface number density in apparent magnitude 
bins. The goal of the present study is to use this informa- 
tion to constrain the structure of the Galaxy parameterised 
in the manner explained above. 

A Monte Carlo method was adopted to simulate the ob- 
servational statistics. Random points were generated based 
on several quantities, including the luminosity function, dis- 
tances weighted by the input density profile, and a Gaussian- 
type observational error, so that the observed star counts 
can be reproduced. For these simulations, we adopted the 
three-dimensional extinction model of the Milky Way de- 
rived from COBE observations (Drimmel, Cabrera-Lavers & 



Lopez-Corredoira 2003). Instead of correcting the observa- 
tional data for the very complicated interstellar reddening, 
we applied reddening corrections to the stellar distributions 
predicted by the model. The extinction. Ay, is in general 
small throughout the areas selected in this study, with a 
maximum of about 0.104 mag in the direction (0°, 55°) and 
a minimum of 0.0385 mag in the anti-Galactic Centre direc- 
tion {b — 80°), at a distance of about 600 pc along the line 
of sight. Av is obtained through linear interpolation for dis- 
tances loss than 600 pc. For greater distances. Ay is taken to 
be constant. The g-to-F-band extinction ratios were taken 
from Girardi et al. (2004), and thus we can also obtain the 
apparent g and r-band magnitudes at each theoretical point. 

Wc can now caluclatc the theoretical projected surface 
number density, theoretical star counts as a function of ap- 
parent g and r magnitude, and theoretical star counts as a 
function of {g — r) colour down to the SDSS limiting mag- 
nitudes, and hence we are now ready to compare these with 
the observational data. 



4 RESULTS AND ANALYSIS 

4.1 Analysis of the first group of observational 

data at 6=60° 

The discrete small sky areas of group 1 , covering a complete 
circle at the Galactic latitude of 6=60° , form a unique probe 
of the structure of the Galaxy in the current SDSS data 
release. Direct measurements of the stellar number density 
distribution within this data set can tell us whether or not 
it is axisymmetryical. 

We tested this idea by fitting the surface density dis- 
tribution along the circle with an axisymmetric halo. Fig. 4 
shows the projected surface number density distributions of 
the theoretical model and the observational distribution for 
stars of g,r £ [15,22] mag, where the diamonds represent 
the observational projected stellar surface number density. 
To account for observational fluctuations, each sky field was 
divided into four subfields of 1° x 1° (except for the area 
[200°, 60°], which is too small for this procedure), so that 
the statistical fluctuations can be measured. The resulting 
fluctuations over the average of the four subfields are shown 
in the flgure as error bars on the data points. The dashed line 
is the theoretical prediction for an axisymmetric model. The 
dashed theoretical line was normalised by taking the aver- 
age values of the surface density at (0°,60°) and (180°, 60°) 
from the observational data. It is worthwhile to point out 
that the two parts of the dashed theoretical line, i.e., that 
between 1 — 0° and 180°, and that between 180° and 360°, 
are not perfectly symmetrical because of extinction. In addi- 
tion, due to the non-uniform extinction, the surface number 
density curves shown in Fig. 4 are not smooth. It is evident 
that the axisymmetric model does not fit the observations 
very well. Comparing the observational data with the ax- 
isymmetric model prediction, it is obvious that the obser- 
vational data between i = 0° and 180° on the one hand, 
and between i = 180° and 360° on the other are highly 
asymmetrical; the projected surface density in the fields at 
i > 180° is systematically higher than that in the fields 
mirror-symmetrically located on the other side of the merid- 
ian. The axisymmetric theoretical value is a little higher 
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than derived from the observations for longitudes from 0° 
to 180°, and much lower than implied by the observational 
data for £ = 180° to 360° . The deviations from the axisym- 
metric model are rather smooth; they seem to represent a 
systematic shift from the theoretical model. 

To find out what kind of stars are responsible for these 
asymmetric statistics, we performed a number of tests. In 
Table 2 we list the star counts of stars with g,r € [15, 22] 
mag in the paired sky fields selected symmetrically with re- 
spect to the 0° meridian plane, as well as the relative fluctu- 
ations (defined as [surface density of £2 - surface density of 
£1]/ [surface density of £1], where ii and I2 are the Galactic 
longitudes of the paired fields). In Table 3, faint stars with 
g,r £ [19,21] mag are considered in the same way. They 
show the same trend as the previous data set, resulting in 
an uneven distribution with an even higher asymmetric rar 
tio. 

In Fig. 3, star counts in four colours, (u — g), [g — r), 
(r — i) and (i — z), in the fields towards the directions of 
(90°, 60°) and (270°, 60°) are compared. The magnitude 
ranges are, respectively, 15 < w ^ 22, 15 ^ g < 22, 15 ^ 
r ^ 22, 15 < « < 21, and 15 ^ ^; ^ 20. In each of the 
colour distribution profiles there are two peaks, as shown in 
Fig. 3 (the blue tail at (w — <?) ~ 0.1 is mostly due to white 
dwarf stars). The left-hand peak is dominated by halo stars, 
while the one on the right is dominated by thin disk stars, as 
inferred from the colours. This is because most of halo stars 
are far more distant than the bulk of the disk stars, and 
therefore only the luminous stars in the halo (predominantly 
main-sequence stars) can be detected. For staxs on the main 
sequence, the intrinsically brighter ones have bluer colors. In 
addition, the halo stars are generally more metal poor, and 
hence tend to populate the blue peak. Based on a similar 
argument related to the stellar metallicity, thin disk stars 
would form the red peak. Therefore, it is clear that the excess 
of the stellar number density in the field at (270°, 60°) with 
respect to that at (90°, 60°) is due to the halo component. 
The amplitude of the excess depends on colour, as shown in 
Fig. 3; it is larger in {u — g),{g — r) and (r — i) than in (i — z). 
It may be caused by the diflferent limiting magnitudes, or it 
might be an indication of the specific properties of the stars 
that causes the excess. 



4.2 Fitting the projected surface number 
densities with a triaxial halo model 

A triaxial halo model, as described in Section 3, is adopted 

to describe the asymmetry found in our data set. Given the 
same limiting magnitudes as those of the observations, the 
theoretical surface densities are calculated using a model of 
which the structural parameters are constrained by a com- 
parison with the observations. 

The three-component model (including the thin and 
thick disks, and the halo) is used to calculate projected sur- 
face number densities in the sky areas corresponding to those 
of the observations. As the main objective of the current 
work is to derive the structure of the halo, for the parame- 
ters of the thin and thick disk we adopt the results of Chen 
et al. (2001). The parameters of the halo are adjusted in 
such a way that the observational data are best fitted. The 
integration of the observational surface density over all of 



the selected fields is used to constrain the surface density of 
the model. 

There are a total of six parameters in our halo model, 
namely the power-law index n, the two axis ratios p and q, 
and the three coordinate rotation angles 5, (j), 9 (based on 
which the x,y,z axes of the triaxial ellipsoidal model can 
be translated to that of the Galaxy). The possible ranges 
of the parameters are estimated using initial tests over a 
preliminary grid, which is given in Table 4. 

Chi-squarcd minimisation tests comparing the theoret- 
ical and observational data sets were carried out using the 
algorithm of Press et al. (1992). For a non-Gaussian distri- 
bution of discrete data points, Pearson's can be expressed 
as, 

where Ri is the theoretical frequency of data points in the i**" 
bin, and 5*^ is the equivalent frequency of the observations. 
N is the number of bins, m the number of model parameters, 
and therefore {N — m) is the number of degrees of freedom. 
Assuming that the errors in the theoretical and observa- 
tional frequencies are distributed in a Poissonian fashion, 
we obtain cr^ = (^/E^ + ^/Sf)^. By changing the model 
parameters, the parameter combination that generates the 
minimum averaged of the 12 groups can be obtained. 

Models for all possible combiuatious of the six param- 
eters were calculated, and the corresponding x^s were eval- 
uated to find the parameter combination resulting in the 
minimum x^ value. It is obvious that the models with (p, q) 
in the range of (0.5, 0.5), (0.6, 0.5), (0.7, 0.6), and (0.8,0.6) 
lead to better values for x^i regardless of any fine-tuning 
of the other four parameters. The models provide only weaJj 
constraints for the power-law index n. When n is fine-tuned, 
the axial ratio should be adjusted accordingly, in order to 
obtain the optimum x^- The overall trend is that for larger 
n, larger p or q ratios are needed. 

Figs 4 and 5 show the observational surface densities in 
the selected regions and those of our theoretical models cal- 
culated by using a triaxial halo with the following parameter 
combination, n = 2.2, p = 0.5, g = 0.5,9 = 60°,^ = -10°, 
and <p = —10°, and the corresponding x^ is 2.02. (This is the 
parameter set that leads to the minimum Xg-n see Section 
4.3 for more details). Fig. 4 shows the results for the group- 1 
fields, spanning a circle at & = 60°. A comparison between 
the solid and dashed lines in Fig. 4 shows that the triaxial 
halo model (the solid line) fits the observations much better 
than an axisymmetric halo (the dashed line) , with the triax- 
ial halo clearly showing the degree of asymmetry as revealed 
by the observations. The only sizable deviations of the tri- 
axial model from the observations are found near £ = 260° 
(where the theoretical prediction for the surface density is 
slightly higher), and near £ — 330° (where the theoretical 
value is somewhat lower). The six panels in Fig. 5, labelled 
A-F, show the surface densities of both the observations (di- 
amonds and triangles) and the models (dotted and dashed 
lines) for groups 2 to 12. Panel A shows the surface densities 
in the fields at £—0° and ^=180° as a function of Galactic 
latitude. Panels B through E show the surface densities in 
the sky-field pairs mirror-symmetrically selected on either 
side of the degree Galactic meridian plane. Panel F de- 
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picts that in the £=150° field, the companion field of which 
is not shown because the data has not yet been released. 
Prom these figures, it is clear that there are larger surface 
densities in the sky areas between i = 180° and 360° than 
in their corresponding mirror-symmetric fields. Those fields 
thus show the same trend as the sky areas in the first group, 
at 6 = 60° . The theoretical model fits the observations quite 
well in general, with several exceptions. For the three groups 
at £=0°,240°, and 270°, the values are higher than those of 
the observations around a latitude of 6=55°. For the group 
at £ = 300° and £=270°, the theoretical values are lower 
than the observational value around 6 = 70° . For the group 
at £ = 330° and £—240°, the theoretical values are lower 
than the observational values around b = 75° and b = 55°, 
respectively. 

In this paper, a power-law model is used to interpret the 
observations for the purpose of clarity. Models implemented 

with more complex physics, such as the Hernquist model 
(Newberg & Yanny 2005b), can fit the observations better. 
This could be one of the reasons that our model does not 
fit the observations perfectly. Modifications to this problem 
will be considered in a future study. 

4.3 Fitting star counts in the {g — r) colour 

In the previous section, the projected surface number den- 
sities were fitted using a triaxial halo model. Star counts 
combined with colour information can, however, reveal more 
information about the stellar populations in the data and 
constrain the model better. Therefore, we performed model 
fits in the [g — r) colour. The model described in Section 3 
is now used to calculate star counts in the (g — r) colour. 

By using Eq. (1), the projected stellar number density in 
each absolute magnitude bin can be obtained. Monte Carlo 
simulations were carried out for absolute magnitude and dis- 
tance, so that artificial stars could be generated for a given 
model, and hence theoretical apparent magnitudes could be 
generated. The {g — r) colour as well as its statistical uncer- 
tainty at each point can then be obtained from the colour- 
magnitude relation. 

Comparisons between the theoretical and observational 
star counts for all the fields were also done following Eq. (4). 
Subsequently, the xfg-r) values of all fields in the 12 groups 
were averaged to obtain the mean value X(j_r)' Intensive 
and elaborate computational tests were done for all possible 
cases. Each set of paramet ers repr esents a model which can 
be evaluated by its Bind xfg-r)- '^''^^ found that and 
X^g-r) exhibit similar trends. However, when X(g_r) reaches 
its minimum value, x^ may not be at its minimum. There- 
fore, we have to sacrifice minimizing the latter in order to 
keep X(g_r) small. Table 5 lists the best-fit parameter com- 
binations for power-law indices n from 2 to 2.6. Because of 
the six-dimensional parameter space, multiple solutions were 
derived. Fitting actual star counts in {g — r) is much more 
difficult than fitting the projected surface number density. 
The minimum value is about 1.735, but that of xfg-r) ^ 
a few times higher. In our model, the power-law index n still 
leaves much room for fine-tuning, from 2.0 to 2.6 (beyond 
which and xfg-r) become too large). If constrained only 
by the projected number density (see Section 4.2), no firm 
constraint on the axial ratios {p,q) can be reached. When 



supplemented with star counts in (g — r), each n has only one 
best-fit {p,q) combination. For n ^ 2.4, {p,q) is (0.5,0.5). 
This best-fit parameter set is different from that of the tra- 
ditional axisymmetric model, even ii p = q. In this model, 
the y axis is equal to the z axis, and therefore resembles 
the rotationally symmetric model, but with an inclination 
and tilt with respect to the Galactic plane. We would like to 
call all these p = q cases also triaxial models. In such cases, 
good fits to the observational data can also be achieved. 
However, the ellipticity seems too oblate. When n — 2.4, 
{p,q) becomes (0.6,0.5); when n = 2.5 or 2.6, (p, g) goes to 
(0.7,0.6). The larger n becomes, the higher the axial ratio 
will need to be. Directly following from the power-law distri- 
bution, the number density decreases more rapidly for larger 
n when going away from the Galactic Centre. A higher axial 
ratio tends to enhance the number density in the directions 
of the axes. Hence, for larger n, a larger axial ratio is re- 
quired to balance the effect of the larger n in order to fit the 
observational data. Such a correlation is the same as that 
for the direct surface number counts in Section 4.2. 

Figs 6 and 7 show star counts in {g — r) for the 6=60° 
fields. The theoretical values were calculated with the same 
model as in Figs 4 and 5. Our purpose is to compare 
the mirror-symmetric pairs; four pairs of fields at approx- 
imately this latitude are considered. Because the sky field 
at (30° , 60° ) is not included in the survey, we chose the field 
at (30°, 65°) to be paired with that at (330°, 65°). The thin- 
line histogram represents the observational data, the solid 
black histogram describes the theoretical results. In Fig. 6, 
panel A shows the results for the sky field at (0° , 60° ) , and 
panel B for (180°, 60°). The panel pairs C-D, E-F, G-H and 
1-3 are for four pairs of mirror-symmetric sky fields. Panel K 
in Fig. 7 is for the field at (150° , 60°), a single field without a 
companion on the other side of the Galactic meridian plane. 
Comparing the observational data in the paired sky fields, 
the histograms for £ from 180° to 360° exhibit higher halo 
star peaks, and similar peaks for the thin disk with respect 
to those on the other side of the meridian. The theoretical 
model follows the same trend as the observations. 

Apart from the problem of the structural parameters, 
the stars in the Galaxy are simplified as three populations 
(thin and thick disks, and halo), and it is further assumed 
that each population internally has the same age and metal- 
licity. All of these assumptions will surely affect the distri- 
bution of the stars in the {g — r) colour. A more precise 
g-band luminosity function and better partitioning of the 
populations are needed to interpret the continuously updat- 
ing large-scale survey data. 



5 DISCUSSION AND CONCLUSION 

Star counts based on the observational data of the Milky 
Way analysed in the present paper show an asymmetry in 
the large-scale areas of the sky surrounding the northern 
Galactic pole, and therefore they refiect a certain non ax- 
isymmetric structure. These observations may be explained 
in one of two ways. The present smooth Galactic halo it- 
self may have an asymmetric structure, which can be repre- 
sented by a triaxial ellipsoid, possibly tracing the dark halo 
gravitational potential that might be triaxial (Jing & Suto 
2002). The triaxiality proposed as such is not unique; it is 
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considered common for dark haloes in extragalactic systems. 
Mazzei & Curir (2001) analysed the effects of a triaxial dark 
halo on the bar-triggered star formation and feedback pro- 
cesses after the formation of the disk. They indicated that 
the dark halo has important effects on the evolution of bary- 
onic matter. The star-formation rate is not only linked to the 
total mass of the dark halo, but also to the dynamical state 
of the dark halo. Is it then normal to form an asymmet- 
ric stellar halo in such an asymmetric gravitational field? In 
the history of formation and cvohition of the Galaxy, the 
halo experienced accretion and merger processes. Wyse & 
Gilmore (2005) claim that merger events and accretion of 
satellite galaxies in the past history of the Galaxy do play 
some role in the evolution of the Milky Way. The outer halo, 
in particular, might bo doiniuatcd by substructures that are 
likely the remnants of interactions. If the halo were very 
sparse originally and were gradually built up from accretion 
and merger events, it would be expected that there have 
been accretion events happening throughout the lifetime of 
the Milky Way. This would allow for any sizable features 
from early times to be smoothed out via dynamical evo- 
lution, leading to the formation of a smooth structure at 
the present time. Presumably, the star streams should have 
a random distribution due to the random infall angle and 
intensity, and therefore - on average - the overall features 
would have a smooth distribution after relaxation. Newberg 
and Yanny (2005b) selected F-type turn-off stellar data in 
eight SDSS stripes from the DR3 to map the distribution 
of stars in real space. They selected fields that most likely 
contain only halo stars. Sky areas near any previously iden- 
tified overdensities caused by star streams were also avoided 
in their studies, but their results can be altered by the Virgo 
overdcrisity discovered later by Juric et al. (2005). After they 
compared the observational data with the results of five tri- 
axial models, they obtained the best-fit parameter set for 
an oblate ellipsoid with its major axis pointed at 50-70° 
from the line of sight towards the Galactic Centre from the 
Sun, and with a 4-6° inclination with respect to the Galac- 
tic plane. The minor axis, z, has a length equivalent to 65% 
of the major axis. The intermediate axis, y, is about 75% of 
the major axis (Newberg & Yanny 2005b). 

An alternative cause of an asymmetric halo may be 
found in large-scale star streams that are embedded in the 
smooth background and that have not yet been melted into 
the halo completely. Also using DR4 data, Juric et al. (2005) 
analysed the structure of the Galaxy using the photometric 
parallax method. They found that there are significant over- 
densities in the range (5 < ^;/kpc < 15), with the overden- 
sity peaking in the direction of Virgo, and the entire overden- 
sity region covers about one thousand square degrees on tlie 
sky. They selected both a sky field in the overdensity region 
and a control sky field mirror-symmetrically located on the 
other side of the meridian with respect to the target field, 
as shown in figure 24 of Juric et al. (2005). By comparing 
stars with 0.2 < {g — r) < 0.8 mag and 18 < r < 21.5 mag 
in both fields, they estimated the surface brightness and lu- 
minosity of the overdensity in the target field as J2r ~ "^2-5 
mag arcsec"^, and Lr = 0.09 x 10^ L©. Based on this result, 
we can calculate the additional number of stars per square 
degree due to the overdensity. The magnitude range used 
by Juric et al. (2005) spanned from r = 18 to 21.5 mag. 
Adopting the distance to the overdensity as 10 kpc, as also 



assumed by Juric et al. (2005), the lower and upper limit 
of the r-band absolute magnitude of their sample stars is 
therefore from 6.5 to 3 mag. From the transformation rela- 
tion r = V - 0.49(S -V) + 0.11 (Chen et al. 2001), the 
limiting magnitudes in the F-band can be obtained (taking 
the V - {B — V) relation for main-sequence stars from Cox 
1999). The Mass - V magnitude relation from Reid (2002), 

logM = 0.477-0.135Mv-hl.228xl0"^xM^-6.734xl0"^My 

(5) 

is used to calculate the mass of the most luminous (1.45 
Mq) and the faintest stars (0.81 M©) corresponding to the 
respective observed magnitude thresholds. Integrating the 
luminosity- mass ratio, weighed by mx IMF (Reid 2002), 
the additional stellar density number per square degree can 
be calculated, 

N = Ltot/(c j (L/m)m(/>(m)dm), (6) 

where c is the normalisation constant of the IMF. The de- 
rived number is 614 stars per square degree, which is indeed 
enough to match the number density excess with respect 
to a symmetric structure. However, this is still not enough 
to fully resolve the problem. By subtracting the CMD of a 
control field from that of the sky area characterised by the 
density excess, the flux of the Virgo star stream was derived. 
However, if a triaxiality does exist in the structure, there 
would be no proper criteria to distinguish the contribution 
due to the star stream from that of triaxiality. 

As demonstrated in Fig. 4, there is a large region of 
stellar underdensity with respect to an axisymmetric model, 
which is centred at about £ = 150° in the direction of Ursa 
Major (UMa) with an amplitude less than the Virgo over- 
density. This was also shown in figure 3 of Newberg & Yanny 
(2005a), where they quote a minimum at £ = 155° for F-type 
turn-off stars, very close to what we have found in this pa- 
per. The size of the UMa underdensity region is smaller than, 
but comparable to, that of the Virgo overdensity. The UMa 
underdensity is also clearly visible from figure 24 of Juric 
et al. (2005; the blue colour-coded region near £ = 150°), 
and more prominently shown on the colour map in a recent 
news release on the SDSS web page. Such an underdensity 
region can also be found in the spatial distribution of M- 
giants from 2MASS observations (figure 27 of Juric et al. 
2005). The Virgo overdensity and the UMa underdensity re- 
gions are well aligned. If the overdensity is to be accounted 
for only by a star stream, what is then the cause for the ob- 
served underdensity aligned with it? This infers that there is 
some substructure within the halo that cannot be accounted 
for only by a star stream such as that of the Virgo over- 
density. As shown in Fig. 4, the overdensity is observed in 
a very large region, from £ = 180° to 360°, which is more 
than a quarter of sky at the latitude of 6 = 60° . It cannot be 
fully understood by a stream covering only 1000 deg^. The 
triaxial halo model may provide a plausible explanation, be- 
cause the distribution of the stellar surface number density 
at 6 = 60°, including these over- and under-densities, can 
be fitted simultaneously by a triaxial halo model. 

To disentangle the effects of a triaxial halo and the 
large-scale star stream, kinematic and chemical information 
is needed. 
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The substructures detected in the SDSS data are sum- 
marised in Table 2 of Newberg, Yanny & Rockosi (2002); two 
of their star streams pass through the sky fields selected in 
this study, namely S341+57-22.5 and 8297+63-20.0. Since 
S341+57— 22.5 is below the magnitude limits of the survey, 
it can hardly alter our results. S297+63— 20.0 may have a 
minor influence on our results because of its limited size in 
the given direction. However, it will not change the general 
large-scale picture of the current work. 

As the primary goal of the current work, a study of the 
large-scale structure of the Galactic halo using star counts 
is presented. The star streams can not easily be corrected 
for from the SDSS data. Therefore, we chose all-sky cover- 
age to try to fit the average structure. We have shown that 
the stellar projected surface density of the halo stars at £ 
from 180° to 360° is significantly higher than that at £ from 
0° to 180° at a latitude of 60° below the SDSS magnitude 
limit. The observed stellar projected surface density distri- 
butions in the sky fields selected for the current work can 
be well fitted by using a triaxial halo model. The parameter 
combinations of the best-fit triaxial models are, power-law 
index n from 2.0 to 2.6, 9 from 55° to 65°, <j>, £, less than 5°, 
{p,q) either (0.5,0.5), (0.6,0.5), or (0.7,0.6). These param- 
eter sets agree with the one of Newberg & Yanny (2005b), 
except that their axial ratios are larger than ours. The re- 
sults of the current work can still be affected by lumps in 
the star streams that have not yet been smoothly dispersed 
into the halo. 

We are left with a number of questions, such as: what 
are the bulk of the stars in the star streams? What is the 
proportion of stars in the streams across the whole halo? 
Do star streams contribute significantly to the formation of 
a triaxial halo, or do they simply act as minor disturbances 
to the overall profile? How many of the star streams have 
already dispersed into the halo and become a smooth back- 
ground, and how many of them can still be picked up from 
the smooth halo background? There is no way to have any 
clear answer for these questions from the present data. Work 
in this subject will leap forward with the coming on-line of 
future projects aimed at spectroscopic sky surveys such as 
SEGUE, LAMOST, SDSS H, GAIA, and further photomet- 
ric surveys of the southern sky. 
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Figure 1. A sketch of the selected sky areas of the current study in Galactic coordinates. 
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Figure 2. A sketch of the orientation of the coordinate frames. 



Table 1. The Galactic coordinates of the selected sky fields in each group 



group (l,b) 



1 {0°,60°), (10° ,60°), (40° ,60°), (50° ,60°), (60° ,60°), (70° ,60°), 
(80° ,60°), (90°, 60°), (100° ,60°), (110° ,60°), (120° ,60°), (130°, 60°), 
(140° ,60°), (150° ,60°), (160° ,60°), (170° ,60°), (180° ,60°), (190° ,60°), 
(200° ,60°), (230° ,60°), (240° ,60°), (250° ,60°), (260°, 60°), (270°,60°), 
(280° ,60°), (290° ,60°), (300° ,60°), (310° ,60°), (320° ,60°), (330° ,60°), 
(340° ,60°), (350° ,60°) 

2 (0°,55°), (0°,60°), (0°,65°), (0°,70°) 

3 (30° ,65°), (30° ,70°), (30°,75°) 

4 (60° ,55°), (60° ,60°), (60° ,65°), (60° ,70°), (60° ,75°), (60° ,80°), (60° ,85°) 

5 (90° ,55°), (90° ,60°), (90° ,65°), (90°, 70°), (90°, 75°), (90° ,80°), (90°, 85°) 

6 (120° ,55°), (120° ,60°), (120° ,65°), (120°, 70°), (120°,75°),(120°,80°),(120°,85°) 

7 (150° ,55°), (150° ,60°), (150° ,65°), (150°,70°),(150°,75°), (150°,80°), (150°,85°) 

8 (180° ,55°), (180° ,60°), (180° ,65°), (180°, 70°) 

9 (240° ,55°), (240° ,60°), (240° ,65°) 

10 (270° ,55°), (270° ,60°), (270° ,65°), (270°, 70°), (270°,75°) 

11 (300° ,60°), (300° ,65°), (300°, 70°), (300°, 75°) 

12 (330°, 60°), (330°, 65°), (330°, 70°), (330°, 75°) 
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Figure 3. Star counts in four colors, namely (it — g), (g — r), (r — i) and (i — z). The thick solid histogram represents the star counts at 
(270°, 60°), and the thin histogram describes star counts at (90°, 60°); the two sky fields were chosen in mirror-symmetry with respect 
to the of 1=0° Gala<;tic meridian plane. 
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Figure 4. Stellar surface number density distribution at a latitude of 60°. The diaiiioiids describe the surface number density of the 
observational data with g and r-band magnitudes from 15 to 22 mag; the error bars represent density fiuctations. The dashed line 
describes the theoretical result calculated based on an axisymmetric halo model, and the solid line describes the theoretical result for 
triaxial halo model. 

Table 2. The relative deviations of the sky field pairs for p,r S [15, 22] mag. (see text for details). 



£l (degree) 


£2 (degree) 


asymmetry ratio(%) 


10.000 


350.000 


7.49773 


40.000 


320.000 


21.3195 


50.000 


310.000 


19.9415 


60.000 


300.000 


11.2650 


70.000 


290.000 


20.6879 


80.000 


280.000 


23.9994 


90.000 


270.000 


18.9822 


100.000 


260.000 


14.3713 


110.000 


250.000 


9.98441 


120.000 


240.000 


16.5926 


130.000 


230.000 


14.8426 


160.000 


200.000 


13.5624 


170.000 


190.000 


5.12844 
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Figure 5. Fits of the surface iiuiiiber density using groups 2-12. The solid and dashed lines are the theoretical predictions, while the 
diamonds and triangles show the observational data. Panel A; solid line, diamond, 1=0°; dashed line, triangle, 1=180°. Panel B: solid 
line, diamond, 1=30°; dashed line, triangle, 1=330°. Panel C: solid line, diamond, 1=60°; dashed line, triangle, 1=300°. Panel D: solid 
line, diamond, 1=120°; dashed line, triangle, 1=240°. Panel E: solid line, diamond, 1=90°; dashed line, triangle, 1=270°. Panel F: solid 
line, diamond, 1=150°. 
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Figure 6. Fits of star count statistics in {g — r) colour, where g and r are both within [22, 15] mag. The thin lines describe the 
observational data, the solid black lines show the theoretical results. Panel A: the results for (0°,60°), Panel B: (180°, 60°), Panel C: 
(30°, 65°) Panel D: (330°, 65°) Panel E: (60°, 60°) Panel F: (300°, 60°). 



Table 3. The relative deviations of the sky field pairs, similar to Table 1^ but for g, r £ [19, 21] mag 



£i (degree) 


£2 (degree) 


asymmetry ratio(%) 


10.000 


350.000 


10.2431 


40.000 


320.000 


32.7824 


50.000 


310.000 


26.8430 


60.000 


300.000 


24.9343 


70.000 


290.000 


35.5734 


80.000 


280.000 


33.3278 


90.000 


270.000 


29.9915 


100.000 


260.000 


26.4650 


110.000 


250.000 


15.3487 


120.000 


240.000 


16.3516 


130.000 


230.000 


12.8257 


160.000 


200.000 


19.2948 


170.000 


190.000 


5.18586 
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Figure 7. -continued. Panel G: (90°, 60°), Panel H: (270°, 60°), Panel I: (120°, 60°) Panel J:(240°, 60°), Panel K: (150°, 60°) 
Table 4. Input ranges of the parameters. 



parameter lower limit upper limit step 



2 
0.5 
0.5 

55° 
-20° 



3 0.1 

0.9 0.1 

0.9 0.1 

65° 5° 

5° 5° 
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Table 5. The best-fit pcirameters for fixed power-law index n at [2,2.6] 



power-law index 


(degree) 


P 


g 


f (degree) 


cp (degree) 




A(g_T.) 


2. 


65. 


0.5 


0.5 


-15. 


-5. 


2.302 


7.092 


2.1 


65. 


0.5 


0.5 


-20. 


-10. 


1.735 


6.842 


2.2 


60. 


0.5 


0.5 


-10. 


-10. 


2.02 


6.632 


2.3 


55. 


0.5 


0.5 


5. 


-15. 


2.869 


6.909 


2.4 


65. 


0.6 


0.5 


5. 


-20. 


2.965 


7.073 


2.5 


65. 


0.7 


0.6 


-15. 


-10. 


3.055 


7.204 


2.6 


65. 


0.7 


0.6 


-15. 


-10. 


3.539 


7.257 



